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^ ' Influence of surrounding matter on the properties of clusters is considered by an approach combin- 

Qing the methods of statistical and quantum mechanics. A cluster is treated as a bound A'^ -particle 
system and surrounding matter as thermostat. It is shown that, despite arbitrary strong interactions 
^Nl ' between particles, cluster energy can be calculated by using the controlled perturbation theory. The 

accuracy of the latter is found to be much higher than that of the quasiclassical approximation. 
Spectral distribution is obtained by minimizing conditional entropy. Increasing the thermostat tem- 
perature leads to the depletion of bound states. The characteristic temperature when bound states 
flj ■ become essentially depleated defines the temperature of cluster evaporation. The inverse process of 

lowering the thermostate temperature, yielding the filling of bound states, corresponds to cluster 
condensation. 
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^ . I. INTRODUCTION 

C^ 

The formation of clusters is a common tendency of matter. Clusters of quarks form nucleons, clusters of nucleons 

form nuclei, clusters of nuclei form atoms and molecules, clusters of the latter form various objects on the Earth, 

'^ . including alive beings. The Earth itself, as well as other planets, are clusters. Planets clusterize around the Sun or 

Ch ■ around other stars. Galaxes in the Universe are also a kind of clusters. 

J^ , Generally, a cluster is a compact group of objects bound together for sufficiently long time. More exactly, the 

lifetime of a cluster is to be much longer than the characteristic times describing its motion as a whole. In the ideal 

case, when the lifetime of a cluster is infinite, the latter is absolutely stable. Usually, clusters have finite lifetimes, 

thus, being metastable, although the lifetime can be extremely long. 

J^ ' In what follows, for concreteness, we shall mean under a cluster a group of TV particles forming a bound state 

IliHi ■ because of attracting forces between them. There are two main approaches of treating clusters theoretically, that of 

f^ ' statistical and of quantum mechanics. The statistical approach permits to define the relative concentrations of different 

rvq , clusters and other thermodynamic functions of a statistical systems, provided that the basic characteristics of each 

T-H type of clusters, such as mass, binding energy, and quantum numbers, like spin, are given a priori. This approach 

C ■ ', describes, for example, the relative abundances of chemical elements in the universe [1]. In the droplet models of 

^^ ■ condensation [1,2] the noninteracting droplets are characterized by their volumes and surface tension. Alternatively, 

the interaction potentials between droplets are to be postulated [3]. The appearance of droplets is typical of the 

liquid-vapour phase transition occurring either in usual gases or in nuclear matter [4-6] . 

Recently, a great attention has been attracted to the use of the statistical approach for describing the deconfinement 

phase transition in quantum chromodynamics. When increasing temperature or density, hadron matter is expected 

r-| ■ to decompose into quark-gluon plasma [7-12]. High density of nucleons can be reached in ultrarelativistic nuclear 

Q , collisions [13-16] and both high density and temperature, in heavy-ion collisions [17-20]. Numerical simulations with 

O lattice models [21-24] undoubtfuUy show that nonperturbative phenomena dominate in a wide range around the 

deconfinement phase transition, up to temperatures of about 2Tc [25,26]. These nonperturbative phenomena are 

thought to be caused by the formation of quark and gluon clusters [27-32] inside quark-gluon plasma, similar to the 

appearance of liquid droplets in fog. The deconfinement-confinement phase transition, to our mind, is an analog of 

the evaporation-condensation phase transition. 

In the quantum-mechanical approach, one considers the properties of an individual cluster in the vacuum and the 

interactions of a small number of clusters. For instance, in last years the properties of multiquark clusters have been 

intensively studied [33-38]. Such clusters, as various scattering experiments indicate [39-43], can exist in nuclei and 

even in deuteron [44,45]. More information on multiquark clusters can be found in reviews [46,47]. Other examples of 

clusters whose description is mainly based on quantum mechanics are nuclei themselves, alpha clusters [48], resonant 

nuclear clusters [49,50] and metallic clusters [51]. 

The properties of a cluster, defined when the latter is in the vacuum, can be strongly changed if the cluster is inserted 

into matter. Moreover, the stability of a cluster can be essentially influenced by surrounding medium. Therefore, 
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when one calculates the concentrations of clusters in a statistical system, fixing cluster characteristics a priori, while 
the latter can be noticeably disturbed by the influence of matter, would yield wrong conclusions. At the same time, 
for understanding how matter influences cluster properties, one needs to have a method of their accurate definition. 
To describe the properties of a cluster is usually not very easy since one has to deal with bound states of strongly 
interacting particles. 

In the present paper we consider the following mutually related questions: (i) How to accurately calculate the 
energy spectrum of bound states for an TV -particle cluster? (ii) What are the statistical weights of such states? (iii) 
How the evaporation or condensation of a cluster is caused by surrounding matter? 

The questions addressed above are, of course, too complicated to be exhaustively answered in one paper. Here, we 
suggest some principal ways of solving the problem and illustrate the ideas by a simplified model: A cluster consists of 
N particles interacting with each other through attractive pair forces of power law; surrounding matter is replaced 
by thermostat. We hope that the suggested approach can be applied to more realistic problems. But first we have to 
be convinced that there exists the meth 

od permitting to answer the above formulated questions in principle. 

II. CLUSTER ENERGY 

Consider a system of N particles defined by the Hamiltonian 

Hn= I i}jH7)K{7)tlj{7)d7 + 

+ - / V (^)V' (^ )^{\ r - r \)'4}{r )ip{r)d r d r , (1) 

in which K{r) is a kinetic-energy operator; tp{''~) is a field operator being a column [V's(^)] m the space of 
internal degrees of freedom, such as spin, isospin and so on. The interaction potential does not depend on internal 
degrees of freedom. 

The average energy of the system can be exactly expressed through the second-order density matrix. There are 
several forms of such expressions [52-54]. In the second-quantization representation, which we use here, the average 
energy can be written as 

1 /■ ^ ^/ ^ ^/ ^ ^ 

{Hn) = 2(N-1) J ^^''^ ^ ^i)'^(^2 - ^2)^2(^1, ?-2)x 

X ^2(7-1, ^2, r^, r^)d r^dr^dr^d r^, (2) 

where the reduced Hamiltonian is 

i^2(^i, r 2) = K{7^) + K{72) + iN- 1)$(| 7i - ^2 I) (3) 

and the second-order density matrix 

P2{7i, 7 2, ^1, ^2) = {i'\~r ^)'4)\7 {)'4){7 i)i}){7 2)) . (4) 

The form (2) can be easily obtained using the fact that, when the number of particles is fixed, the second- and 
first-order density matrices are connected by the relation 

P2(7^1, ^2, ^1, 72)d "^2= {N - l)pi(7i, T^J, 

pi(r, r ) = (■0T(r )-0(r)). 

Let us note that there is another way, not mentioned in literature, of deriving (2). Write the Hamiltonian in the 
Schrodinger representation 

N N 

H{7,, 72,...,7^)^Y1 ^(^i) + 2 E *(l ^^ - ^j- !)■ (5) 

i=l i^j 



Passing to the Heisenberg representation, one transforms the first and second terms in (5) to the first and second 
terms in (1), respectively. However, nothing prohibits us to rewrite (5) in the identical form 

1 " ^ ^ 

in which H2 is given by (3). From the latter expression, in the Heisenberg representation we obtain 

Hn = 2(N - X\ I vi^Wir )H2{r, r )^p{r )tp{r)d r d r . (6) 

And from (6) one immediately gets (2). 

The second-order density matrix (4) can be presented as an expansion 

P2iri, r2, r^, r ^) ^ 2^S„^*a(ri, r2)-^*^{r^, r ^) (7) 

a/3 

over an orthonormalized basis {^q} • The expansion coefficients, because of the normalization condition 

obey the sum rule 

J2Bc.c. = N{N-1). (8) 

a 

With the notation 

the average energy (2) becomes 

It is natural to choose the basis {^q} in the expansion (7) as a set of eigenfunctions of the reduced Hamiltonian 
(3), 

i/2(r 1, 72)*„(7i, ^2) = s„*„(7i, 72). (10) 

Equation (10) describes quantum mechanics of a pair of particles. Such a pair can be called [54] a pairon. Defining 



Pa 



N{N - 1)' 



Ep" = 1' (11) 



we see that Pa has the meaning of the probability that a pairon is in a state a . Now, the average energy (9) reduces 
to 

{Hn) ^ -Yp-^^o.. (12) 

a 

For concreteness, let us take the kinetic-energy operator in the standard form 

V2 



K{r) = -, 



2mo ' 
TTig being a particle mass. Then, using the center-of-mass and relative coordinates 

-^ 7i + ^2 ^ ^ ^ 

i?= ?; , ^ = ^1-^2, 



and introducing the masses 

M=2mo, fi^—, 

one can, as is known, factorize the eigenfunction in (10) into the product 

in which ^g(_R) is a plane wave, with a wave vector Q , corresponding to a pairon center-of-mass motion; Yim is 
a spherical harmonic; and the pairon relative motion is described by a function Xni satisfying the equation 

HN{r)xni{^) = Eni{N)xni{r) (13) 

with the relative Hamiltonian 

The pairon energy E^ in (10) reads 



*2 



Ec^^ + E„iiN). (15) 

Thus, the pairon quantum index a is specified as a set a = {Q, 71,1,771} of six indices: three components in the 

wave vector Q= {Qa \a = 1, 2, 3} , each Qa G R ; the radial quantum number n , which is n = 0, 1, 2, . . . for bound 
states and n € R for unbound states; the orbital quantum number I = 0, 1,2, . . . and the azimuthal quantum 
number m ~ —I, —l + l,..., —1, 0,1, ... ,1 — 1,1 . 

Now we need to define the statistical weights of the energy levels (15). We assume that the considered system of 
N particles is immersed into thermostat with temperature T . To find the weights Pa in (12) we can proceed as 
follows. Define the entropy 

S = -^Pa\npa. 

a 

The weights Pa have to satisfy the normalization condition in (11) and give an average energy per particle according 
to (12), that is 



^Pa = l, 1^(^n) = ^PaEa. 



N 

a a 

In addition, there can exist some constraints or selection rules which usually can be related to operators / being 
integrals of motion whose averages can be written as 



{I) =^PaIa- 



For instance, / can be a projector to a space of a given symmetry. Introduce the notation 

TTa =exp(-7/a), 

where 7 is a Lagrange multiplier. Minimizing the conditional entropy 



S^ S + P (^{Hn) -J2po^eA + 
il-lnZ)lj2p^-l\ +7f(/)-5IP"^")' (16) 



with the Lagrange multipliers (3 ^ T ^ , 7 , and 1 — InZ , we get 

Pa^ -^exp{~/3Ea); Z = ^7r„ exp (-/3Sa) . (17) 

a 

As an example of selection rules defining tTq consider the case when the second-order density matrix possesses the 
symmetry property 

P2(^2, ri,r-^,r^) :^ ±p2{r i, r2,r^,r^). 

Respectively, the functions $„ in the expansion (7) are either symmetric or antisymmetric with respect to coordinate 
permutation, 

*a(»^2, »^l) =±*a(^l, ?^2). 

The permutation ri^^r2 is equivalent to the transformation 

R^R, r^-r, ■d^ir-'d, ip^ip + n. 

Using the property of spherical functions 

we have 

*a(?^2,^l) = (-l)'*a(^l, ?^2). 

Therefore the coefficients Bap in (7) should have the form 



BaP ~ TT^TTijCaP, 



in which 



+ ^ +_ , 1; / = 0,2,4... 
" ~ ' "1 0; i = 1,3,5,... 



0; /==0,2,4... 
^" -^' ^ ^ 1; ? = 1,3,5,... 



In this case, tt^ is a projector to even or odd states, respectively. Then, from (17) we derive 

V-- ±2? + i r /? , ,] 

Pnl = 2^PQnlm = TT; — CXp <^ --Enl{N) ^ ; 



Z,„t=^^±(2/ + l)exp|-^^„,(iV)|. (18) 

nl ^ ^ 



In this way, the average energy (12) takes the form 



{Hn)^^NT + E{N), (19) 



where 

EdV] = 



E{N) = -Y,PniEni{N). (20) 



nl 



The interpretation of the obtained result is straightforward. Since we are thinking of the considered system as of an 
N -particle cluster, we may say that the first term in (19) is the kinetic energy of the cluster, and (20) is the internal 
energy of particles inside the cluster. A part of the spectrum Eni{N) must correspond to bound states. In order 
that we could, in principle, speak about an N -particle cluster, the number of bound states must be not less than 

N . 



III. BOUND STATES 

To be able to form bound states, the particles should interact with each other through forces containing an attractive 
part. The interaction potential in the reduced Hamiltonian (3) includes the factor (N — 1) . When N ^ 1 , then, 
even if the potential <& is small, the product {N — !)<!> can be large. Moreover, the interaction between particles 
forming bound states is practically never small. All this means that the standard perturbation theory will not work 
for describing bound states. However, the problem is treatable by the controlled perturbation theory [55-58] which is 
applicable to arbitrary strong interactions. 

For demonstrational purpose, let us take an attractive potential of power law 

$(r) = -Ar-", (21) 

in which A, ly > and r e R_|_ . Then the Hamiltonian (14) reads 

1 ^^ 1(1 + ^) ,^r -,^ ^ 
2/i dr^ 2jir^ 

The latter, using a scaling, can be reduced to the dimensionlcss form 



HN{r) = --^ + -\^ -{N-1) -. (22) 



Id' m + i}_9 ,,„. 

The return from (23) to (22) is made by substituting 

to 
r ^ y/JJUJr, g—^N-1. (24) 

As zero approximation we may accept the Hamiltonian 

Id' l{l + l) u 
^'--2d^ + ^^-^ (^^) 

with the Coulomb interaction including yet unknown parameter u . The eigcnfunctions of (25) are 



xi"Ar) 



{n + 2l + iy. 



1/2 , / o \ ' + 1 

' 1 / 2ur ^ 



X exp ; L 



n+l+1 \n+l+l 
2ur 



X 



-2J+1 



n + l + 1 J " \n + l + lj ' 

where Ll^ is an associate Laguerre polynomial. The eigenvalues of (25) are 

u' 
Eo{n,l,9,u)^^ ^^^^^^^^^ . (26) 

Starting from (25), we can invoke the Rayleigh-Schrodinger perturbation theory with the perturbation 

In doing this, we encounter the integral 

JL(-) ^ f\''+'+'^e-'L'J+\t)Lf+\t)dt. (27) 

The latter can be presented as a sum. For instance, if s — n , then 

,,, (n + 2/ + l)!,^(-l)JT(p + 2/ + 3 + i/) , ^ , 

^">) = nl ]2 plin-p)Wip + 2l + 2f -^-P-'-''^^ 



where 



In particular, 



Vn[x) = , , ; ipo[x) = 1. 
r{x) 



jLi-l) = ^'' + ^\+^^\ jU--) = m + 3 - 



For the first approximation we have 



ti oil 

Ei{n,l,g,u) ^Eo{n,l,g,u) + - — - - — ' Ani{i^); (28) 

[n + I + I)'' (n + t + 1)^ 



_ nUU-^) ( 2 -"-^ 

Anl — 



here 

(n + 2? + l)! V"- + ^ + l 
_ r(2^ + 'i-v) ( 2 ' ""^ 

^°' - (2; + !)! Iv^TT 

The procedure can be continued yielding E]^{n,l,g,u) . Then, at each step k , the unknown parameter u is to be 
transformed to a function Uk(n,l,g) , such that the sequence {e-k{n,l,g)}'^^Q of terms 

efc(n, /, g) = Ek{n, I, g, Uk{n, I, g)) (29) 

would be convergent. Because of their role of forcing the convergence of the approximation sequence, the functions 
Mfc(n, /, g) are called control, or governing, functions. Their explicit form is defined by fixed-point conditions. One of 
the simplest such conditions is the condition of minimal sensitivity [59] , 

-^Ek{n,l,g,uk) = 0. (30) 

OUk 

When equation (30) has no solution for the control function Uk , one can use other types of fixed-point conditions 
[55-58]. If there exist several solutions for Uk , than one has to resort to the stability analysis [57,60] choosing that 
solution which makes the procedure more stable. 

In the first approximation we have only one solution for the control function 

ui{n,l,g)^[vAr,i{v)g]^'^'^-''K (31) 

From (28) and (29) we find 

ei(nj,g) = — - — ^ ■ 32) 

Eq.(32) shows that for v > 2 there are no bound states since all energy levels become positive moving to the 
continuous spectrum. When v = 1 , then 

ui{n,l,g)=g, A„i(l) = 1 (j^ = 1), 

and (32) yields 

^^("'^'"^^' 2(n// + l)^ ('^ = ^)' 

that is, we get an exact spectrum corresponding to the Coulomb potential. 
For asymptotically small powers :/ — > , since 

4„(0) = 2(n + / + i) (" + 2^ + l)! ^ ^^^(0) = {n + l + 1)^, 



from (31) wc have 
Therefore (32) gives 

In the case of n = , (32) becomes 



ui(n, Z,5) ~ (n + Z + l)^/pF (ly —> 0). 



\imei{n,l,g) = -g. 



ei{0, 1,9) = - 



2-u 
8iy 



2iyT{2l + 3-iy) 



ii + i)i2i + iy. 



2/(2-,.) 



And the ground-state energy level, when n = / = , is 

ei(0,0,g) = -^ 



[2!.r(3 - iy)9f^^^''''> . 



(33) 



(34) 



We calculated also the second-order approximation 62(^1, /, g) and the corrected first approximation, following the 
method of refs. [57,58,60]. The latter approximation has the form 



e*i{n,l,g) = ei(n, /,g) exp{Ai(n, 0(5*}, 



(35) 



in which Xi(n,l) is a local Lyapunov exponent for the approximation flow and (5* is a damping parameter. However, 
we shall limit ourselves by the first approximation (32) because the expressions for the second and for the corrected 
first approximations are essentially more cumbersome, althougt, as numerical analysis has shown, they do not differ 
much in accuracy from the first approximation. 

The stability of the calculated procedure in the controlled perturbation theory [55-58,60] can be checked in the 
following way. Define the coupling function gk{n, I, f) by the equation 



Eo{n,l,g,UkinJ,g)) ^ f. 
This, in the case of (26) and (31), leads to the equality 

ui{n,l,g) ^ {n + I + 1)^/^. 
For the coupling function we obtain 



(36) 



9i{nJJ) 



{n + l + lf-'' 



(-2./) 



1-1//2 



Now we have to define the points of the trajectory of the approximation cascade as 

yk{n, I, /) = ek{n, I, gk{n, I, /)). 
The starting point is yoin, /, /) = / and the first point is 



The mapping multipliers are 



so that 



A necessary stability condition is 



yi{n,lj) = /. 



of 



f^i{n,l,.f) = 



2-v 



(37) 



(38) 



(39) 



(40) 



(41) 



\^^k{n,l,f)\ < 1, 



(42) 



which shows the existence of contraction at the k -step with respect to the starting point. As we see, (41) satisfies 
(42) if i^ > 1 . 

Continuing the stabihty analysis for the second point of the cascade trajectory, we come to the mapping multiplier 

^i2{nJ,f)^fii{nJJ) + Xi{n,l). (43) 

Numerical calculations evince that the maximal Lyapunov exponent 

Ai = supAi(n, (44) 

is positive and small, < Ai ^ 1 . Whence, the maximal mapping multiplier (43) is 

2-1/ 



snp\n2inJJ)\ 



+ Ai 



The stability condition (42) for the mapping multiplier (43) is uniformly valid for all n, I if 

.>^-^.l + f. (45) 

Another stability condition requires that the Lyapunov exponents Afc would be negative [58,60]. 
To check the accuracy of the procedure, let us introduce the notation 

ek{n,l) = ek{n,l,l), (46) 

using which (29) can be written as 

ek{nJ,g)^ek{nJ)g'^^^-'\ (47) 

The energy levels (46) can be compared with those obtained by direct numerical solution of the corresponding 
Schrodinger equation [61,62]. 

The Schrodinger equation with the Hamiltonian (23) can also be solved in the quasiclassical approximation [63] 
which yields 

ewKB{n,l,g) 



TT 2(2 - v){n + l + l) + {v-l){2l + 1) 

There exists a generalization of the quasiclassical approximation based on the use of modified Airy functions [64] , but 
this technique, as it seems, is useful only for one-dimensional problems. 
Introducing the notation 

ewKB{n,l) = ewKB{n,l,l), (49) 

we can cast (48) into 

ewKB(n,l,g) = ewKB{n,l)g^/'^^^''\ (50) 

The form (49) is convenient to compare with direct numerical solutions [61,62]. 

The most proper characteristic of the precision of an approximate method is, apparently, the relative error. Thus, 
the error of (46) is defined as a percentage deviation 

_ ek{n,l) -e{n,l) 

Sk{n, I) = — — 100% 

|e(n,0| 

from an exact numerical value e(n, I) . The percentage errors of (35) and (48) are defined in the same way, as 
£^(ri, /) and £wKB{n, I) , respectively. The uniform accuracy of an approximation, with respect to all energy levels, 
is characterized by the maximal error 

Sk = sup \ek{n,l)\. 



Analogously, the maximal errors e^ and Ewkb are defined. 

The accuracy of the controlled perturbation theory and that of the quasiclassical approximation, for power-law 
potentials, is demonstrated in tables I-III, as compared to the exact numerical values e(n, I) and with each other. 
In the table I the results are presented for the potential (21) with the power i> = 0.2 . The maximal errors are 
ei = 0.8% and Swkb = 4% . The errors for the second approximation 62 (n,/) and for the corrected first 
approximation e\{n,l) are not displycd since they are very close to ei(n, Z) . The accuracy of (35) in the case of 
the damping parameter (5* = 1 is practically the same as in the case of i5* = 1/2 , for the latter case being slightly 
better. The accuracy, when passing from ei(n, /) to 62 (n,/) , does not improve because for v = 0.2 the stability 
condition (45) is not valid. The corrected approximation e\{n,l) is not more accurate than ei{n,l) since the 
Lyapunov exponent (44) is positive. 

Table II presents the calculations for v = 0.8 . The maximal errors are £1 = 0.3% and Ewkb — 6% . The 
accuracy of 62 and ej is pretty much the same as of ei . The maximal error e^ practically does not change 
when replacing S^ — 1 by (5, — 1/2 . The absence of the improvement in accuracy is again caused by the fact that 
the stability conditions do not hold for u = 0.8 . Nevertheless, as in the previous case of j^ = 0.2 , the controlled 
perturbation theory is much more accurate than the quasiclassical approximation. 

Table III contains the results for v > 1 . For v = 1 , all approximations become exact answers. For 1/ > 1 the 
stability condition (45) now holds, and the accuracy of 62 is slightly higher than of ei . The Lyapunov exponent 
(44), as earlier, is positive, so e* docs not improve ei . Again, the quasiclassical approximation is essentially less 
precise than the controlled perturbation theory. 

IV. EVAPORATION TEMPERATURE 

For the considered case of the power-law attraction (21), the spectrum of the relative motion of particles inside a 
cluster, defined by eq.(13), has the form 

E„i{N)=Lue{nJ)g^^^'-'\ (51) 

where 

^=(^■^^2)1/(2-^), ^^!^, 9 = N-l. 

The value of e(n, I) , as is shown in the previous section, can be calculated, with a good accuracy, by the controlled 
perturbation theory [55-58]. For example, we can use ei{n,l) from (32). 
Define the ground-state level of (51) as 

E^{N) = E^o{N) {n = l = Q). (52) 

Here we should substitute 

ei(0,0) = -^ [2z.r(3 - j.)]2/(2--) . (53) 

The internal energy of particles inside a cluster is given by (20) with the spectral weights (18). The whole spectrum 
has two parts, discrete and continuous. The discrete spectrum corresponds to bound states with negative energies 
(51). The continuous spectrum is related to unbound states with positive energies. 

When temperature increases from zero, the internal energy of a cluster (20) increases from a negative value to zero 
and then to positive values. The characteristic temperature, Tc , at which the cluster ceases to exist, can be defined 
as the temperature where E{N) becomes zero, 

E{N)^0 (T = T,). 

This temperature may be named the evaporation temperature. The absolute value |_B(Af)| is called the binding 
energy. So, the evaporation temperature is that where the binding energy is zero. 

Approximately, the evaporation temperature can be obtained as follows. Note that when the temperature of 
thermostat is such that 2T <C |ii^o(A^)| ^ then the lower bound states are mainly occupied, while the unbound 
states of continuous spectrum arc practically empty. This situation describes a well defined cluster. Recall that the 
thermostat here models surrounding matter, and temperature can be interpreted as a measure for the kinetic energy 
of collisions. If the latter are intensive, so that 2T ^ |ii^o(A^)| j then the bound states are strongly depleted and the 

10 



unbound states are essentially populated. This means that the cluster itself is actually destroyed, since a cluster is 
well defined only if the particles are bound inside it. The characteristic temperature 

UN)^^\Eo{N)\ (54) 

specifying the desintegration of a cluster can be called the evaporation temperature, or the condensation temperature 
if the inverse process of lowering temperature is considered. The energy Eo(N) in (54), defined in (52), is the 
ground-state energy of a pair of particles from the collection of N particles. This explains the appearance of the 
factor i in (54). 

It is convenient to introduce the dimensionless evaporation temperature 

CO 

For the latter, taking account of (51) and (53), wc have 

*-(^) = ^ [2^r(3 - i^){N - l)]2/(2--) . (56) 

In the case of the purely Coulomb interaction (ly ~ 1) , (56) reduces to the expression 

tciN)^\iN~lf (z. = l). 

Then (54) is recovered by multiplying the latter by lu with 

a^ 2 ^ r T\ 

to — — mnc : a = —- (ly = L). 

if A = e'^ , where e is the electron charge, then a = 1/137 is the fine-structure constant. In the latter case, 
LO ^ 13.613ey , or 15.796 x lO^K . 

The evaporation of a cluster is the manifestation of its instability against excitations caused by the influence of 
surrounding. There can exist another instability limiting the size of a cluster. This is related to the definition of the 
cluster mass 

Mat = moiV + -^ > 57 

as of a nonnegative quantity. Here E(N) is the internal energy of particles inside the cluster, being given by (20). 
Due to the dependence of E{N) on temperature, the mass (57) is also dependent on T , so that Mn increases 
together with T . Let us introduce a function "f{T) by the relation 

N 
E{N) EE -jiT)EoiN), (58) 

in which 

So(iV)=we(0,0)(iV-l)2/(2-.)_ 

At zero temperature, < 7(0) < 1 . With the notation (58) the cluster mass (57) reads 

Mn = moN + 7V^^^^e(0, 0){N - 1)2/(2-'^), (59) 

2c^ 

As far as for T < Tc the function 7(T) is positive, and e(0, 0) is always negative, the second term in (59) is 
negative. Therefore, with increasing N , the cluster mass (59) decreases and becomes zero at A^ = No{T) which is 



No{T) ^ 1 + 
In the case of v — 1 , since e(0, 0) = —1/2 , we have 



2moc2 ^'-"^^ 



^7(T)e(0,0) 



(60) 



^°(^)^i+V^^^^ ^^^'^- 
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If we assume that a = 1/137 and take j{T) « 1/2 , then Nq ss 549 . 

Eq.(60) defines the maximal number of particles in a cluster. If N surpasses A^o , the cluster desintegrates into 
several smaller clusters. With temperature increasing from zero to Tc , the maximal number of particles (60) increases. 
When T^Tc-0 , then j{T) ^ . If we suppose that j{T) ^T-T^ as T ^ T^ , then No{T) - (T - Tji-'^/^ 
That is, the maximal number of particles in a cluster tends to infinity with the critical index 1 — z^/2 . The function 
7(r) plays here the role of an order parameter. 

The evaporation of a cluster, though reminds a phase transition, is, of course, not the same. Generally, we can 
consider a cluster of any number of particles N > 1 . Even we may take N = 2 . It would be too bold to speak about 
a phase transition for two particles. An actual evaporation-condensation phase transition occurs in an ensemble of 
many clusters, or droplets, which interact with surrounding as well as with each other. However, if the number of 
particles forming a cluster is very large, say N ^ oo , then the evaporation of such a cluster can have many features 
of a phase transition. For example, a crystal can be treated as a cluster. The transition of the particles of this 
cluster from bound to unbound states happens at the melting temperature Tm which is an analog of the evaporation 
temperature. The temperature T„i can be estimated as follows [65] . According to the meaning of the transition from 
bound to unbound states, the related characteristic temperature is to be close to the ground-state binding energy, 
Tjn ~ \Eo\ . The latter is about the effective depth of the interaction potential |-Eo| ~ 'I'm ■ So, 

where C is a constant of the order of one. Take, e.g., the elements of the inert group interacting through the 12 — 6 
Lennard-Jones potential. Then $m = e , where e is a constant of the Lennard- Jones potential. Estimates give 
C = 1/V2 . The formula 

T - — 

"~ V2 

evaluates the melting temperatures surprisingly well, which is demonstrated in table IV for Ne, Ar, Kr and Xe 
at zero pressure: The calculated Tm is very close to the experimental value T!^'^ . 

V. CONCLUSION 

A cluster is defined as an object by N particles in bound states. For any, arbitrary strong, interactions between 
particles the spectrum of bound states can be calculated with a good accuracy by means of the controlled perturbation 
theory. The latter is about an order more accurate than the quasiclassical approximation. 

The properties of a cluster in the vaccuum and inside matter can be essentially different. The differences are traced 
by analysing a simple model when surrounding matter is replaced by thermostat. With increasing temperature, the 
cluster begins to evaporate. At the characteristic temperature of evaporation, Tc , the cluster desintegrates into its 
constituents. The evaporation temperature is approximately equal to one half of the ground-state binding energy. 

The number of particles in a cluster is limited by iVo depending on the interaction potential between particles and 
the thermostat temperature. The limit iVo increases with temperature, tending to infinity as T ^^ Tc ■ 

The term evaporation is applicable to clusters of any number of particles A^ > 1 . The evaporation of small clusters 
is rather a continuous process than an abrupt transformation. The evaporation of large clusters can acquire features 
of phase transitions. 
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Table Captions 

Table I. The accuracy of the energy levels calculated using the controlled perturbation theory and the quasiclassical 
approximation, for the interaction potential with v — 0.2 . 

Table II. The same as in table I, but for the potential with the power v ~ 0.8 . 

Table III. The accuracy of the ground-state energy found by using the controlled perturbation theory and the 
quasiclassical approximation, for potentials with different powers. 

Table IV. The estimation of the melting temperature at zero pressure, for the inert-group crystals with the 
Lennard-Jones interaction. 
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Table I 



n 


I 


e{n,l) 


ei(%) 


ewKBm 


£wkb/£i 






1 
2 


-0.72526 
-0.63319 
-0.58215 


0.50 
0.29 
0.19 


-1.2 

-2.8 
-3.5 


2.4 
9.7 
18 


1 



1 
2 


-0.60834 
-0.56730 
-0.53733 


-0.32 
0.19 
0.24 


-0.40 
-1.3 
-2.0 


1.3 

6.8 
8.3 


2 



1 
2 


-0.55191 
-0.52680 
-0.50628 


-0.76 
-0.04 
0.14 


-0.22 
-0.83 
-1.4 


0.29 
21 
10 



Table II 



n 


I 


e(n,0 


ei(%) 


£wkb{%) 


\£wkb/£i\ 






1 
2 


-0.48336 
-0.19843 
-0.11707 


0.21 
0.05 
0.20 


-2.0 
-4.9 

-5.8 


9.5 
98 
29 


1 



1 
2 


-0.18334 
-0.11151 
-0.07739 


-0.10 
0.29 
0.21 


-0.58 
-2.2 
-3.5 


5.8 
7.6 
17 


2 



1 
2 


-0.10517 
-0.07421 
-0.05635 


-0.27 
-0.08 
0.17 


-0.18 
-1.5 
-2.3 


0.67 
19 
14 



Table III 



V 


e(0,0) 


ei(%) 


e2(%) 


S-WKB^y^ 


\£WKBle2\ 


1 

1.25 
1.5 


-0.5 
-0.69938 
-2.3687 




1.4 
12 




1.2 
11 




7.1 

28 


1 

5.9 
2.5 



Table IV 





e{K) 


Tm = ^K) 


t:^^p{k) 


TVe 


36 


25 


24 


Ar 


121 


86 


84 


Kr 


163 


115 


117 


Xe 


232 


163 


161 
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